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Abstract 

Abstract. We discuss geometric properties of non-Noether sym- 
metries and their possible applications in integrable Hamiltonian sys- 
tems. Correspondence between non-Noether symmetries and conser- 
vation laws is revisited. It is shown that in regular Hamiltonian sys- 
tems such a symmetries canonically lead to a Lax pairs on the alge- 
bra of linear operators on cotangent bundle over the phase space. 
Relationship between the non-Noether symmetries and other wide 
spread geometric methods of generating conservation laws such as bi- 
Hamiltonian formalism, bidifferential calculi and Frolicher-Nijenhuis 
geometry is considered. It is proved that the integrals of motion as- 
sociated with the continuous non-Noether symmetry are in involu- 
tion whenever the generator of the symmetry satisfies a certain Yang- 
Baxter type equation. Action of one-parameter group of symmetry 
on algebra of integrals of motion is studied and involutivity of group 
orbits is discussed. Hidden non-Noether symmetries of Toda chain, 
nonlinear Schrodinger equation, Korteweg-de Vries equations, Benney 
system, nonlinear water wave equations and Broer-Kaup system are 
revealed and discussed. 
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1 Introduction 



Symmetries play essential role in dynamical systems, because they usually 
simplify analysis of evolution equations and often provide quite elegant solu- 
tion of problems that otherwise would be difficult to handle. In Lagrangian 
and Hamiltonian dynamical systems special role is played by Noether sym- 
metries — important class of symmetries that leave action invariant and 
have some exceptional features. In particular, Noether symmetries deserved 
special attention due to celebrated Noether's theorem, that established cor- 
respondence between symmetries, that leave action functional invariant, and 
conservation laws of Euler-Lagrange equations. This correspondence can be 
extended to Hamiltonian systems where it becomes tighter and more evi- 
dent then in Lagrangian case and gives rise to Lie algebra homomorphism 
between Lie algebra of Noether symmetries and algebra of conservation laws 
(that form Lie algebra under Poisson bracket). 

Role of symmetries that are not of Noether type was suppressed for quite a 
long time. However after some publications of Hojman, Harleston, Lutzky 
and others (see [16], [36], [39], [40], [49]- [54]) it became clear that non-Noether 
symmetries also can play important role in Lagrangian and Hamiltonian dy- 
namics. In particular according to Lutzky [51], in Lagrangian dynamics there 
is definite correspondence between non-Noether symmetries and conservation 
laws. Moreover, unlike Noether's case, each generator of non-Noether sym- 
metry may produce whole family of conservation laws (maximal number of 
conservation laws that can be associated with non-Noether symmetry via 
Lutzky's theorem is equal to the dimension of configuration space of La- 
grangian system). This fact makes non-Noether symmetries especially valu- 
able in infinite dimensional dynamical systems, where potentially one can 
recover infinite sequence of conservation laws knowing single generator of 
non-Noether symmetry. 

Existence of correspondence between non-Noether symmetries and conserved 
quantities, raised many questions concerning relationship among this type 
of symmetries and other geometric structures emerging in theory of inte- 
grable models. In particular one could notice suspicious similarity between 
the method of constructing conservation laws from generator of non-Noether 
symmetry and the way conserved quantities are produced in either Lax the- 
ory, bi-Hamiltonian formalism, bicomplex approach or Lenard scheme. It 
also raised natural question, whether set of conservation laws associated with 
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non-Noether symmetry is involutive or not, and since it appeared that in 
general it may not be involutive, the need of involutivity criteria, similar 
to Yang-Baxter equation used in Lax theory or compatibility condition in 
bi-Hamiltonian formalism and bicomplex approach, emerged. It was also 
unclear how to construct conservation laws in case of infinite dimensional 
dynamical systems where volume forms used in Lutzky's construction are no 
longer well defined. Some of these questions were addressed in papers [11]- 
[14], while in the present review we would like to summarize all these issues 
and to provide some samples of integrable models that possess non-Noether 
symmetries. 

Review is organized as follows. In first section we briefly recall some as- 
pects of geometric formulation of Hamiltonian dynamics. Further, in second 
section, correspondence between non-Noether symmetries and integrals of 
motion in regular Hamiltonian systems is discussed. Lutzky's theorem is 
reformulated in terms of bivector fields and alternative derivation of con- 
served quantities suitable for computations in infinite dimensional Hamilto- 
nian dynamical systems is suggested. Non-Noether symmetries of two and 
three particle Toda chains are used to illustrate general theory. In the subse- 
quent section geometric formulation of Hojman's theorem [36] is revisited and 
some samples are provided. Section 4 reveals correspondence between non- 
Noether symmetries and Lax pairs. It is shown that non-Noether symmetry 
canonically gives rise to a Lax pair of certain type. Lax pair is explicitly 
constructed in terms of Poisson bivector field and generator of symmetry. 
Sample of Toda chains are discussed. Next section deals with integrability 
issues. Analog of Yang-Baxter equation that, being satisfied by generator 
of symmetry, ensures involutivity of set of conservation laws produced by 
this symmetry, is introduced. Relationship between non-Noether symme- 
tries and bi-Hamiltonian systems is considered in section 6. It is proved 
that under certain conditions, non-Noether symmetry endows phase space of 
regular Hamiltonian system with bi-Hamiltonian structure. We also discuss 
conditions under which non-Noether symmetry can be "recovered" from bi- 
Hamiltonian structure. Theory is illustrated by sample of Toda chains. Next 
section is devoted to bicomplexes and their relationship with non-Noether 
symmetries. Special kind of deformation of De Rham complex induced by 
symmetry is constructed in terms of Poisson bivector field and generator 
of symmetry. Samples of two and three particle Toda chain are discussed. 
Section 8 deals with Frolicher-Nijenhuis recursion operators. It is shown 
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that under certain condition non-Noether symmetry gives rise to invariant 
Frolicher-Nijenhuis operator on tangent bundle over phase space. The last 
section of theoretical part contains some remarks on action of one-parameter 
group of symmetry on algebra of integrals of motion. Special attention is 
devoted to involutivity of group orbits. 

Subsequent sections of present review provide samples of integrable models 
that possess interesting non-Noether symmetries. In particular section 10 re- 
veals non-Noether symmetry of n-particle Toda chain. Bi-Hamiltonian struc- 
ture, conservation laws, bicomplex, Lax pair and Frolicher-Nijenhuis recur- 
sion operator of Toda hierarchy are constructed using this symmetry. Further 
we focus on infinite dimensional integrable Hamiltonian systems emerging in 
mathematical physics. In section 11 case of nonlinear Schrodinger equation 
is discussed. Symmetry of this equation is identified and used in construction 
of involutive infinite sequence of conservation laws and bi-Hamiltonian struc- 
ture of nonlinear Schrodinger hierarchy. Section 12 deals with Korteweg-de 
Vries and modified Korteweg-de Vries equations. Non-Noether symmetries 
of these equations produce infinite number of conserved quantities in invo- 
lution. The same symmetries give rise to bi-Hamiltonian structure of KdV 
hierarchies. Next section is devoted to non-Noether symmetries of integrable 
systems of nonlinear water wave equations, such as dispersive water wave 
system, Broer-Kaup system and dispersiveless long wave system. Last sec- 
tion focuses on Benney system and its non-Noether symmetry, that appears 
to be local, gives rise to infinite sequence of conserved densities of Benney 
hierarchy and endows it with bi-Hamiltonian structure. 



2 Regular Hamiltonian systems 

The basic concept in geometric formulation of Hamiltonian dynamics is no- 
tion of symplectic manifold. Such a manifold plays the role of the phase 
space of the dynamical system and therefore many properties of the dynami- 
cal system can be quite effectively investigated in the framework of symplectic 
geometry. Before we consider symmetries of the Hamiltonian dynamical sys- 
tems, let us briefly recall some basic notions from symplectic geometry. 
The symplectic manifold is a pair (M, u) where M is smooth even dimen- 
sional manifold and w is a closed 
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duo = 



(1) 



and nondegenerate 2-form on M. Being nondegenerate means that contrac- 
tion of arbitrary non-zero vector field with uo does not vanish 

i xU j = X = (2) 

(here i x denotes contraction of the vector field X with differential form). 
Otherwise one can say that uo is nondegenerate if its n-th outer power does 
not vanish (oo n ^ 0) anywhere on M. In Hamiltonian dynamics M is usu- 
ally phase space of classical dynamical system with finite numbers of degrees 
of freedom and the symplectic form oj is basic object that defines Poisson 
bracket structure, algebra of Hamiltonian vector fields and the form of Hamil- 
ton's equations. 

The symplectic form uo naturally defines isomorphism between vector fields 
and differential 1-forms on M (in other words tangent bundle TM of symplec- 
tic manifold can be quite naturally identified with cotangent bundle T*M). 
The isomorphic map ^ from TM into T*M is obtained by taking contrac- 
tion of the vector field with uo 

^:X^ -i x u (3) 

(minus sign is the matter of convention). This isomorphism gives rise to 
natural classification of vector fields. Namely, vector field Xh is said to 
be Hamiltonian if its image is exact 1-form or in other words if it satisfies 
Hamilton's equation 

i Xh u + dh = (4) 

for some function h on M. Similarly, vector field X is called locally Hamil- 
tonian if it's image is closed 1-form 

i x uj + u = 0, du = (5) 
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One of the nice features of locally Hamiltonian vector fields, known as Li- 
ouville's theorem, is that these vector fields preserve symplectic form uj. In 
other words Lie derivative of the symplectic form ui along arbitrary locally 
Hamiltonian vector field vanishes 

L x u = <^ ix^> + du = 0, du — (6) 

Indeed, using Cartan's formula that expresses Lie derivative in terms of con- 
traction and exterior derivative 

Lx = ixd + di x (7) 

one gets 

L x uj = i x duj + dixOJ = d%x^> (8) 

(since duj = 0) but according to the definition of locally Hamiltonian vector 
field 

dixw = —du = (9) 

So locally Hamiltonian vector fields preserve uj and vise versa, if vector field 
preserves symplectic form ui then it is locally Hamiltonian. 
Clearly, Hamiltonian vector fields constitute subset of locally Hamiltonian 
ones since every exact 1-form is also closed. Moreover one can notice that 
Hamiltonian vector fields form ideal in algebra of locally Hamiltonian vector 
fields. This fact can be observed as follows. First of all for arbitrary couple 
of locally Hamiltonian vector fields X, Y we have L x u = L Y oo = and 

LxLyuJ — LyLx^ = L[x,Y] u = (10) 

so locally Hamiltonian vector fields form Lie algebra (corresponding Lie 
bracket is ordinary commutator of vector fields). Further it is clear that 
for arbitrary Hamiltonian vector field and locally Hamiltonian one Z one 
has 

L z uo = (11) 
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and 



i Xh oo + dh = (12) 

that implies 

Lz(ix h ^ + dh) = L[ Z ,x h ]UJ + ix h L z u + dL z h = L [z ,x h \u + dL z h = (13) 

thus commutator [Z, X h ] is Hamiltonian vector field X Lzh , or in other words 
Hamiltonian vector fields form ideal in algebra of locally Hamiltonian vector 
fields. 

Isomorphism $ w can be extended to higher order vector fields and differential 
forms by linearity and multiplicativity. Namely, 

$ w (XA7)=$ w (X)A$ w (y) (14) 

Since $ w is isomorphism, the symplectic form uj has unique counter image 
W known as Poisson bivector field. Property du = together with non 
degeneracy implies that bivector field W is also nondegenerate (W n ^ 0) 
and satisfies condition 

[W,W}=0 (15) 

where bracket [, ] known as Schouten bracket or supercommutator, is actu- 
ally graded extension of ordinary commutator of vector fields to the case of 
multivector fields, and can be defined by linearity and derivation property 

[Ci A C 2 A ... A C n , S x A S 2 A ... A S n ] = 
{-l) p+q [C p , S q ] A d A C 2 A ... A C p A ... A C n A S x A S 2 A ... A S q A ... A 5 n (16) 

where over hat denotes omission of corresponding vector field. In terms of 
the bivector field W Liouville's theorem mentioned above can be rewritten 
as follows 

[W(u),W] = 0<£> du = (17) 

for each 1-form u. It follows from graded Jacoby identity satisfied by Schouten 
bracket and property [W, W] = satisfied by Poisson bivector field. 
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Being counter image of symplectic form, W gives rise to map $vkj transform- 
ing differential 1-forms into vector fields, which is inverted to the map $ w 
and is defined by 

<& w : u -> W(u); $>w®uj = id (18) 
Further we will often use these maps. 

In Hamiltonian dynamical systems Poisson bivector field is geometric object 
that underlies definition of Poisson bracket — kind of Lie bracket on algebra 
of smooth real functions on phase space. In terms of bivector field W Poisson 
bracket is defined by 

{/) 9} — W(df A dg) (19) 

The condition [W, W] = satisfied by bivector field ensures that for every 
triple (/, g, h) of smooth functions on the phase space the Jacobi identity 

{/{<?, h}} + {h{f, g}} + {g{h, /}} = 0. (20) 

is satisfied. Interesting property of the Poisson bracket is that map from 
algebra of real smooth functions on phase space into algebra of Hamiltonian 
vector fields, defined by Poisson bivector field 

/ - X f = W(df) (21) 

appears to be homomorphism of Lie algebras. In other words commutator of 
two vector fields associated with two arbitrary functions reproduces vector 
field associated with Poisson bracket of these functions 

[Xf,X g ] = X {Lg} (22) 

This property is consequence of the Liouville theorem and definition of Pois- 
son bracket. Further we also need another useful property of Hamiltonian 
vector fields and Poisson bracket 

{/, g} = W(df A dg) = u(X f A X g ) = L Xf g = -L Xg g (23) 
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it also follows from Liouville theorem and definition of Hamiltonian vector 
fields and Poisson brackets. 

To define dynamics on M one has to specify time evolution of observables 
(smooth functions on M). In Hamiltonian dynamical systems time evolution 
is governed by Hamilton's equation 



d_ 
~dt 



f = {h,f} 



(24) 



where h is some fixed smooth function on the phase space called Hamiltonian. 
In local coordinate frame bivector field W has the form 



d d 
W = W bc — A — 



(25) 



and the Hamilton's equation rewritten in terms of local coordinates takes the 
form 



W bc — 

dz c 



(26) 



Note that functions W ab are not arbitrary, to ensure validity of [W, W] = 
condition they should fulfill restriction 



E 

a=l L 



w dW c , u/ dW b , u/ dW bc 

W ab — + W ac — + W ad 



dz n 



dz n 



dz n 



= 



(27) 



and in the same time determinant of matrix formed by functions W ab should 
not vanish to ensure that Poisson bivector field W is nondegenerate. 



3 Non-Noether symmetries 

Now let us focus on symmetries of Hamilton's equation (24). Generally 
speaking, symmetries play very important role in Hamiltonian dynamics due 
to different reasons. They not only give rise to conservation laws but also 
often provide very effective solutions to problems that otherwise would be 
difficult to solve. Here we consider special class of symmetries of Hamilton's 
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equation called non-Noether symmetries. Such a symmetries appear to be 
closely related to many geometric concepts used in Hamiltonian dynamics 
including bi-Hamiltonian structures, Frolicher-Nijenhuis operators, Lax pairs 
and bicomplexes. 

Before we proceed let us recall that each vector field E on the phase space 
generates the one-parameter continuous group of transformations g a = e aLE 
(here L denotes Lie derivative) that acts on the observables as follows 

9a(f) = e aL -(f) = f + aL E f + \{aL E ff + ... (28) 

Such a group of transformation is called symmetry of Hamilton's equation 
(24) if it commutes with time evolution operator 

!*</)= <») 

in terms of the vector fields this condition means that the generator E of the 
group g a commutes with the vector field W(h) = {h, }, i. e. 



[E, W(h)\ = 0. (30) 

However we would like to consider more general case where E is time de- 
pendent vector field on phase space. In this case (30) should be replaced 
with 



-E=[E,W(h)}. (31) 

Further one should distinguish between groups of symmetry transformations 
generated by Hamiltonian, locally Hamiltonian and non-Hamiltonian vector 
fields. First kind of symmetries are known as Noether symmetries and are 
widely used in Hamiltonian dynamics due to their tight connection with con- 
servation laws. Second group of symmetries is less interesting but locally they 
also lead to conservation laws. While third group of symmetries that further 
will be referred as non-Noether symmetries seems to play important role in 
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integrability issues due to their remarkable relationship with bi-Hamiltonian 
structures and Frolicher-Nijenhuis operators. Thus if in addition to (30) the 
vector field E does not preserve Poisson bivector field [E, W] ^ then g a is 
called non-Noether symmetry. 

Now let us focus on non-Noether symmetries. We would like to show that the 
presence of such a symmetry essentially enriches the geometry of the phase 
space and under the certain conditions can ensure integrability of the dynam- 
ical system. Before we proceed let us recall that the non-Noether symmetry 
leads to a number of integrals of motion. More precisely the relationship 
between non-Noether symmetries and the conservation laws is described by 
the following theorem. This theorem was proposed by Lutzky in [51]. Here 
it is reformulated in terms of Poisson bivector field. 

Theorem 1. Let (M, h) be regular Hamiltonian system on the 2n-dimensional 
Poisson manifold M. Then, if the vector field E generates non-Noether sym- 
metry, the functions 



W k A W n ~ k 

y(k) = fc=l,2,...n (32) 

W n 

where W = [E, W], are integrals of motion. 
Proof: By the definition 

W k a W n ~ k = Y {k) W n . (33) 

(definition is correct since the space of 2n degree multivector fields on In 
degree manifold is one dimensional). Let us take time derivative of this 
expression along the vector field W(h), 

^-W k AW n ' k = (^Y^)W n + Y^[W(h),W n ] (34) 

or 



k(^-W) A W k - X A W n ~ k + (n- k)[W(h), W] A W k A W n ^ k ^ 1 = 

(JjL 

(^Y {k) )W n + nY {k) [W(h), W\ A W n ~ l (35) 

(JjL 
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but according to the Liouville theorem the Hamiltonian vector field preserves 
W i. e. 

^-W=[W(h),W] = (36) 
hence, by taking into account that 

l tE= ^ E+mh) , E] = (37) 



we get 



jW = j t [E, W] = [j t E, W] + [E[W(h), W}} = 0. (38) 



and as a result (35) yields 



d Y W W n = (39) 
dt v 1 



but since the dynamical system is regular (W n ^ 0) we obtain that the 
functions are integrals of motion. 

Remark 1. Instead of conserved quantities Y^ l \..Y^ n \ the solutions c\...c n 
of the secular equation 

(W - cW f = (40) 

can be associated with the generator of symmetry. By expanding expression 
(40) it is easy to verify that the conservation laws Y^ can be expressed in 
terms of the integrals of motion C\...c n in the following way 

Y^ = ^J^ £ c n c n ...c lk (41) 

h<t2<---<ik 

Note also that conservation laws Y^ can be also defined by means of sym- 
plectic form u using the following formula 
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while c\...c n conservation laws can be derived from the secular equation 



(L E co - ao) n = (43) 

However all these expressions fail in case of infinite dimensional Hamiltonian 
systems where the volume form 

n = u n (44) 

does not exist since n = oo. But fortunately in these case one can define 
conservation laws using alternative formula 



C« = i wk {L E u) k (45) 

as far as it involves only finite degree differential forms (L E uj) k and well de- 
fined multivector fields W k , further we will use this expression in construc- 
tion of infinite sequence of conservation laws in Korteweg-De Vries, modified 
Korteweg-De Vries and nonlinear Schrodinger equations. Note that in finite 
dimensional case the sequence of conservation laws is related to families 
of conservation laws and in the following way 



C {k) = E c n c l2 ...c lk = . n ;„ ; , r (fc) (46) 
n<ik..<i k {n-k)\k\ 

Note also that by taking Lie derivative of known conservation along the 
generator of symmetry E one can construct new conservation laws 

jY = L Xh Y = =► jL E Y = L Xh L E Y = L E L Xh Y = (47) 
since [E, X h ] = 0. 

Remark 2. Besides continuous non-Noether symmetries generated by non- 
Hamiltonian vector fields one may encounter discrete non-Noether symme- 
tries — noncannonical transformations that doesn't necessarily form group 
but commute with evolution operator 

J t 9(l)=4 t S) (48) 
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Such a symmetries give rise to the same conservation laws 



q(W) k A W n ~ k 

y (fe) = 9{vvj a fc = i > 2,...n (49) 

Sample. Let M be R 4 with coordinates Z\,Zv,Zz,z± and Poisson bivector 
field 



Iir 3 9 d d . n . 

w =^ A ^ + ^ A aJI (50) 

and let's take 



h = \z\ + \z\ + e z ^ (51) 

This is so called two particle non periodic Toda model. One can check that 
the vector field 

E =T, E *£- (52) 



a=l dz a 



with components 



Ex = \z\ - e 23 " 24 -Uz! + z 2 )e Z3 ~ Zi 

£ 2 = \z\ + 2e 23 " 24 + Uz x + 2 2 )e 23 " 24 

E, = 2 Zl + Iz 2 + i(^ + e 23 - 24 ) 

^4 = ^2-^1 + ^(^ 2 2 + e 23 " 24 ) (53) 

satisfies (31) condition and as a result generates symmetry of the dynamical 
system. The symmetry appears to be non-Noether with Schouten bracket 
[E, W] equal to 
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B B B B 

W = [E, W] = zi^— A ^— + z 2 -^— A + 

Bz x Bz 3 Bz 2 BZ4 

B B B B 

oz 1 Bz 2 Bz 3 BZ4 

calculating volume vector fields W k A W n ~ k gives rise to 

„, TTT n d B B B 
WAW = -2— A-A-A- 

oz\ Bz 2 Bz 3 BZ4 

rJr rrr / \ 9 B B B 

W AW — -{ Zl + z 2 — A — — A — — A — — 

Bz\ Bz 2 Bz 3 BZ4 

r\ r\ r\ r\ 

W A W = -2( Zl z 2 - e zs ^ 4 )^ A ^— A ^— A ^— (55) 
and the conservation laws associated with this symmetry are just 



WAW 1, 

It is remarkable that the same symmetry is also present in higher dimensions. 
For example in case when M is R 6 with coordinates 

z l, z 2, Z3, Z4, Z$, Zq (57) 

Poisson bivector equal to 



9 B B B B B . oN 

w = a7 1 ^ i + e^ A a^ + a7 3 ^ (58) 

and the following Hamiltonian 



h = ]-zl + -z 2 2 + \z\ + e Zi ~ Zh + e* 5 " 26 (59) 

jL Zi Zi 
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we still can construct symmetry similar to (53). More precisely the vector 
field 

E = T. E ^ (60) 

a=l uz a 

with components specified as follows 



Ex = l ^z\ - 2e 24 " 25 - ^(zi + z 2 )e Z4 ~ Zh 
E 2 = ]-z\ + 3e 24 " 25 - e 25-26 + Uz l + z 2 )e Z4 ~ Z5 

E 3 = l -z\ + 2e 25 " 26 + l -{z 2 + z 3 )e 25 " 26 (61) 



^4 = 3zi + ^ 2 + \z 3 + |(*? + e 24 ~ 25 ) 
E 5 = 2z 2 - l - Zl + iz 3 + ^ 2 2 + e 24 " 25 + e 25 " 26 ) 

E 6 = z 3 - \z x - X -z 2 + U^zl + e 25 " 26 ) (62) 

satisfies (31) condition and generates non-Noether symmetry of the dynam- 
ical system (three particle non periodic Toda chain). Calculating Schouten 
bracket [E, W] gives rise to expression 



-i d d d d d d 

W = [E, W] = Zl — A — + z 2 — A — + z 3 — A — + 

OZ 1 OZ4 OZ 2 OZ 5 OZ 3 OZq 

d d d d 

C z *- z -> /\ _|_ c 2 5~26 /\ _|_ 

(9^i dz 2 dz 2 dz 3 

d d d d d d . . 

Volume multivector fields W k AW n ~ k can be calculated in the manner similar 
to R 4 case and give rise to the well known conservation laws of three particle 
Toda chain. 
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Y^- l (z + z Iz) WAWAW 



^(2) !/ za-zk z* 2 m W AW AW 

Yi) = 3^ + Z ^ + ^ 3 " - 6 " 6 ) = H/AH/AH/ 

m # A A # 

y( 3 ) = ^ 3 - z 3 e Z4 "* 5 - zie* 8 "* = — — — — (64) 

W AW AW 



4 Non-Liouville symmetries 

Besides Hamiltonian dynamical systems that admit invariant symplectic form 
uj, there are dynamical systems that either are not Hamiltonian or admit 
Hamiltonian realization but explicit form of symplectic structure uj is un- 
known or too complex. However usually such a dynamical systems possess 
invariant volume form Q which like symplectic form can be effectively used in 
construction of conservation laws. Note that volume form for given manifold 
is arbitrary differential form of maximal degree (equal to the dimension of 
manifold). In case of regular Hamiltonian systems, n-th outer power of the 
symplectic form uj naturally gives rise to the invariant volume form known 
as Liouville form 

Q = uj n (65) 

and sometimes it is easier to work with Q rather then with symplectic form 
itself. In generic Liouville dynamical system time evolution is governed by 
equations of motion 

j t f = X(f) (66) 
where X is some smooth vector field that preserves Liouville volume form Q 

^-Q = L x tt = (67) 
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Symmetry of equations of motion still can be defined by condition 

!*</> = («) 

that in terms of vector fields implies that generator of symmetry E should 
commute with time evolution operator X 

[E,X] = (69) 

Throughout this chapter symmetry will be called non-Liouville if it is not 
conformal symmetry of Q, or in other words if 

L E Vt ^ ctt (70) 

for any constant c. Such a symmetries may be considered as analog of non- 
Noether symmetries defined in Hamiltonian systems and similarly to the 
Hamiltonian case one can try to construct conservation laws by means of 
generator of symmetry E and invariant differential form Q. Namely we have 
the following theorem, which is reformulation of Hoj man's theorem in terms 
of Liouville volume form. 

Theorem 2. Let (M, X, Q) be Liouville dynamical system on the smooth 
manifold M. Then, if the vector field E generates non-Liouville symmetry, 
the function 



is conservation law. 
Proof: By the definition 



L E il = JVl. (72) 



and J is not just constant (again definition is correct since the space of volume 
forms is one dimensional). By taking Lie derivative of this expression along 
the vector field X that defines time evolution we get 

LxLeQ = L{x,e]Q + LeLxQ = 

L x (Jty = (L x J)tt + JL x tt (73) 
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but since Liouville volume form is invariant Lx& = and vector field E is 
generator of symmetry satisfying [E, X] = commutation relation we obtain 



(L x J)tt = 



(74) 



or 



Tt J = L * J 



(75) 



Remark 3. In fact theorem is valid for larger class of symmetries. Namely 
one can consider symmetries with time dependent generators. Note however 
that in this case condition [E, X] = should be replaced by 

^E=[E,X\ (76) 

Note also that by calculating Lie derivative of conservation law J along gen- 
erator of the symmetry E one can recover additional conservation laws 

jM = {L E ) m Q (77) 



Sample. Let us consider symmetry of three particle non periodic Toda chain. 
This dynamical system with equations of motion defined by the vector field 



X = 



,24 — 2:5 



a 



24—25 



d d 

^2 5 — 2 6 \ |_ gZ5-26 

dz 2 dz 3 
Odd 

+ Z 1T{ 1" Z 2T{ 1" Z 3T{ — 

OZ4 OZ5 OZq 



possesses invariant volume form 

VL = dz\ A dz 2 A dzs A dz 4 A dz 5 A dze 



(78) 



(79) 



One can check that Q is really invariant volume form, i.e. Lie derivative of 
Q along X vanishes 



d_ 
~dt 



dX 1 dX 2 dX 3 dX 4 dX 5 dX 6 



dzi dz 2 



dz* 



Oza 



8zk 



0Z R 



n = (80) 
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The symmetry (61) is clearly non-Liouville one as far as 



L E Q 



dE 1 dE 2 dE 3 dE A dE b dE 6 



dz\ dz 2 dz 3 dzi dz 5 dz e 
(zi + z 2 + z 3 )dzi A dz 2 A e£z 3 A dz A A dz 5 A dz 6 = 

(zi + z 2 + z 3 )Q 



(81) 



and main conservation law associated with this symmetry via Theorem 2 is 
total momentum 



J = 



z l + z 2 + z 3 



(82) 



Other conservation laws can be recovered by taking Lie derivative of J along 
generator of symmetry E, in particular 

J« = L E J = 

111 

~2 i 2 I 2 I 24—25 1 Z5—ZQ 

2^i " l " 2 2 2 



J( 2 ) = L E JW 
3 3 
4) + 2^ Zl + ^)e 24_25 + -(zi + ^)e" 5 " 
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(83) 



5 Lax Pairs 

Presence of the non-Noether symmetry not only leads to a sequence of con- 
servation laws, but also endows the phase space with a number of interesting 
geometric structures and it appears that such a symmetry is related to many 
important concepts used in theory of dynamical systems. One of the such 
concepts is Lax pair that plays quite important role in construction of com- 
pletely integrable models. Let us recall that Lax pair of Hamiltonian system 
on Poisson manifold M is a pair (L, P) of smooth functions on M with values 
in some Lie algebra g such that the time evolution of L is given by adjoint 
action 
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-L = [L,P] = -ad P L 



(84) 



where [, ] is a Lie bracket on g. It is well known that each Lax pair leads 
to a number of conservation laws. When g is some matrix Lie algebra the 
conservation laws are just traces of powers of L 



j(k) = !L Tr (L k 



(85) 



since trace is invariant under coadjoint action 



d_ 
It 



/(*) 



2d~t 



Tr(L k ) 



-Tr(^L k ) 
2 K dt ' 



^Tr(L k - 1 [L,P])= 1 -Tri[L k ,P]) 







(86) 



It is remarkable that each generator of the non-Noether symmetry canonically 
leads to the Lax pair of a certain type. Such a Lax pairs have definite 
geometric origin, their Lax matrices are formed by coefficients of invariant 
tangent valued 1-form on the phase space. In the local coordinates z a , where 
the bivector field W, symplectic form uj and the generator of the symmetry 
E have the following form 



W 



Wab Q Za A Q Zb 



u ab dz a A dz b E 



ab 



ab 



a 9z a 



(87) 



corresponding Lax pair can be calculated explicitly. Namely we have the 
following theorem: 

Theorem 3. Let (M, h) be regular Hamiltonian system on the 2n-dimensional 
Poisson manifold M. Then, if the vector field E on M generates the non- 
Noether symmetry, the following 2n x 2n matrix valued functions on M 



L a b — X! L ° ad 

dc 



'dW db dE d dE b 

E r — W bc —~ + W dc 



ab 



dz c 



dz c 
dW bc dh 
dz a dz c 



+ W, 



be 



dz c _ 

d 2 h 

dz a z c 



(88) 
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form the Lax pair (84) of the dynamical system (M, h). 

Proof: Let us consider the following operator on a space of 1-forms 

Re(u)=$u([E,<S> w (u)])-L e u (89) 

(here §w and Q w are maps induced by Poisson bivector field and symplectic 
form). It is remarkable that R E appears to be invariant linear operator. First 
of all let us show that Re is really linear, or in other words, that for arbitrary 
1-forms u and v and function foperatorRE has the following properties 

R E (u + v) = R e (u) + R e (v) (90) 

and 

R E (fu) = fR E (u) (91) 

First property is obvious result of linearity of Schouten bracket, Lie derivative 
and maps $vk, 3>w Second property can be checked directly 

R E {fu) = ® U ([E, $ w (fu)}) - L E (fu) = 
$ U ([E, f$w(u)\) - (L E f)u - fL E u = 
$ML E mw(u)) + $ u (f[E, ®w(u)\) - (L E f)u - fL E u = 
L E f^w(u) + f$ u ([E, $w(u)}) - (L E f)u - fL E u = 

$w(u)\) - L E u) = f~R E (u) (92) 

as far as $^$v^( M ) = u. Now let us check that R E is invariant operator 

jTl E = L Xh R E = L Xh {<&uL E §w - L E ) = §u>L [Xh ,E}§w ~ L { x h ,E\ = (93) 

because, being Hamiltonian vector field, X h commutes with maps ®u 
(this is consequence of Liouville theorem) and commutes with E as far as E 
generates the symmetry [X h , E] — 0. In the terms of the local coordinates 
Re has the following form 

— d 

R E = J2 L *bdz a ® t- (94) 

ab ° Z b 
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and the invariance condition 



Jt RE ~ 1 



W(h)RE — 



yields 



d — d ^ T , 9 x 

dt RE= dt^ LabdZ ^dz- b =i: 



ab ab I 

d 

+ L ab (L w{h) dz a ) ® — + Y L ab dz a 
ab ° Zh 

d 

■Kb dz a <g) — + 

abed 



d L 
Jt ab 



d 

dz a ® ^~ 

dz b 



W(h) 



E 

ab 



dt 



JL 

dz b 
d 



d _ r dw ad dh , 

dz a <g> — + L ab — — T—dz c <g> \ 

dZ b abed dZ C OZ d dZ b 



E L abW ad 

abed 



d 2 h , 

a a dz c 
dz c dz d 



d_ 

dz h 



+ 



^ T dW cd dh , d ^ T Txr d 2 h , d 

2^ L ab~^ -K—dZa <g> h 2^ L a bW cd -^—^— dz a ® — 



abed 



E 

ab 



~^L a b + E(^ ac -^ c6 _ LacPcb) 



8 

dz a ® — = 



or in matrix notations 



5 L = ^- 



(95) 



(96) 



(97) 



So, we have proved that the non-Noether symmetry canonically yields a Lax 
pair on the algebra of linear operators on cotangent bundle over the phase 
space. 

Remark 4. The conservation laws (85) associated with the Lax pair (88) 
can be expressed in terms of the integrals of motion c$ in quite simple way: 



I {k) = l -Tr{L k ) = Yc 1 l (98) 
/ i 

This correspondence follows from the equation (40) and the definition of the 
operator Re (89). One can also write down recursion relation that determines 
conservation laws 1^ in terms of conservation laws 
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m—1 

j(m) + (_i)^ m c (m) + ]T = o (99) 

fc=l 

Sample. Let us calculate Lax matrix of two particle Toda chain associated 
with non-Noether symmetry (53). Using (88) it is easy to check that Lax 
matrix has eight nonzero elements 



/ z 1 

z 2 e Z3 ~ Z4 

1 Z! 

V -1 



while matrix P involved in Lax pair 

dt 



[L,P] 



,,23—24 




Z2 



(100) 



(101) 



has the following form 













V e 



23-24 _ p 23-24 



1 OX 

1 



J 



(102) 



The conservation laws associated with this Lax pair are total momentum and 
energy of two particle Toda chain 



/« = -Tr{L) = Zl + z 2 



1. 



1^ = -Tr(L 2 ) = zi + z\ + 2e 23 " 24 
2 



(103) 



Similarly one can construct Lax matrix of three particle Toda chain, it has 
16 nonzero elements 



24 



with matrix P 



Zl 



















z 2 










_gZ 5 -Z6 








Z3 





gZ5-26 








-1 


-1 


Zl 








1 





-1 










1 


1 











^3 



(104) 














1 





\ 
















1 






















1 




pZ4— Z 5 


gZ4-Z5 


















g24~ 25 g z 5— 26 


g25~2 6 











V 





g25~ 2 6 


g25~ 26 








0/ 



(105) 



Corresponding conservation laws reproduce total momentum, energy and 
second Hamiltonian involved in bi-Hamiltonian realization of Toda chain 



JW = l -Tr{L) = Zl + 



Z2 



/( 2 ) = ^Tr{L 2 ) = z\ + z\ + zl + 2e 24 " Z5 + 2e Z5 " 26 



/(3) = -Tr(L 3 ) = z\ + 4 + zl + 3(zi + z 2 )e 24 -* 5 + 3(z 2 + z 3 )e 



25—26 



(106) 



6 Involutivity of conservation laws 

Now let us focus on the integrability issues. We know that n integrals of 
motion are associated with each generator of non-Noether symmetry, in the 
same time we know that, according to the Liouville- Arnold theorem, regu- 
lar Hamiltonian system (M, h) on In dimensional symplectic manifold M is 
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completely integrable (can be solved completely) if it admits n functionally 
independent integrals of motion in involution. One can understand functional 
independence of set of conservation laws Ci,c 2 ...c n as linear independence of 
either differentials of conservation laws dci, dc2---dc n or corresponding Hamil- 
tonian vector fields X C1 , X C2 ...X Cn . Strictly speaking we can say that conser- 
vation laws Ci,C2-..c n are functionally independent if Lesbegue measure of 
the set of points of phase space M where differentials dci, dc 2 ...rfc n become 
linearly dependent is zero. Involutivity of conservation laws means that all 
possible Poisson brackets of these conservation laws vanish pair wise 

{ Ci , Cj } = i,j = l...n (107) 

In terms of the vector fields, existence of involutive family of n functionally 
independent conservation laws Ci,C2--.c n implies that corresponding Hamil- 
tonian vector fields X Cl , X C2 ...X Cn span Lagrangian subspace (isotropic sub- 
space of dimension n) of tangent space (at each point of M). Indeed, due to 
property (23) 

{c i ,c j } = u(X Ci ,X c .) = (108) 

thus space spanned by X C1 , X C2 ...X Cn is isotropic. Dimension of this space is 
n so it is Lagrangian. Note also that distribution X C1 , X C2 ...X Cn is integrable 
since due to (22) 

[X Ct ,X Cj ] = X {ei „ } = (109) 

and according to Frobenius theorem there exists submanifold of M such that 
distribution X Cl , X C2 ...X Cn spans tangent space of this submanifold. Thus 
for phase space geometry existence of complete involutive set of integrals of 
motion implies existence of invariant Lagrangian submanifold. 
Now let us look at conservation laws Y^\Y^ 2 \..Y^ associated with gener- 
ator of non-Noether symmetry. Generally speaking these conservation laws 
might appear to be neither functionally independent nor involutive. How- 
ever it is reasonable to ask the question - what condition should be satisfied 
by the generator of the non-Noether symmetry to ensure the involutivity 
(-[yWj y( m )J. — 0) of conserved quantities? In Lax theory situation is very 
similar — each Lax matrix leads to the set of conservation laws but in general 
this set is not involutive, however in Lax theory there is certain condition 
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known as Classical Yang-Baxter Equation (CYBE) that being satisfied by 
Lax matrix ensures that conservation laws are in involution. Since involutiv- 
ity of the conservation laws is closely related to the integrability, it is essential 
to have some analog of CYBE for the generator of non-Noether symmetry. 
To address this issue we would like to propose the following theorem. 
Theorem 4. If the vector field E on 2n-dimensional Poisson manifold M 
satisfies the condition 



and W bivector field has maximal rank (W n ^ 0) then the functions (32) 
are in involution 



Proof: First of all let us note that the identity (15) satisfied by the Poisson 
bivector field W is responsible for the Liouville theorem 



that follows from the graded Jacoby identity satisfied by Schouten bracket. 
By taking the Lie derivative of the expression (15) we obtain another useful 
identity 



L E [W, W] = [E[W, W]] = [[E, W}W] + [W[E, W}} = 2[W, W] = 0. (113) 



W]]W] = 



(110) 




(111) 



[W, W]=0 & L WU) W = [W(f), W]=0 



(112) 



This identity gives rise to the following relation 



[W,W]=0 & [W(f),W] = -[W,W(f)] 



(114) 



and finally condition (110) ensures third identity 



[W, W}=0 



(115) 



yielding Liouville theorem for W 
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[W,W]=0 [W(f),W]=0 (116) 

Indeed 

[W, W] = [[E, W]W] = [[W, E]W] = 

-[[E, W}W] = —[[E[E, W}}W] = (117) 

Now let us consider two different solutions q ^ Cj of the equation (40). By 
taking the Lie derivative of the equation 

(w - aW) n = o (lis) 

along the vector fields W(cj) and W(cj) and using Liouville theorem for W 
and W bivectors we obtain the following relations 

(W - Cl W) n -\L w{C]) W - {c v c t }W) = 0, (119) 

and 

(w - awy-^CiL^w + {c v Cl }.W) = 0, (120) 

where 

{cj, cj}. = W(dd A dcj) (121) 

is the Poisson bracket calculated by means of the bivector field W. Now 
multiplying (119) by q subtracting (120) and using identity (114) gives rise 
to 

({d, Cj }. - c t {c t , Cj })(W - Cl W) n - x W = (122) 

Thus, either 

{cj, Cj}. - Ci{ci, Cj} = (123) 
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or the volume field (W — CiW) n ~ l W vanishes. In the second case we can 
repeat (119)-(122) procedure for the volume field (W — CiW) n ~ 1 W yielding 
after n iterations W n = that according to our assumption (that the dy- 
namical system is regular) is not true. As a result we arrived at (123) and 
by the simple interchange of indices i <-> j we get 

{cj, cj}, - Cj{ci, cj} = (124) 

Finally by comparing (123) and (124) we obtain that the functions q are in 
involution with respect to the both Poisson structures (since q ^ Cj) 

{q, cj}, = {ci, cj} = (125) 

and according to (41) the same is true for the integrals of motion Y^ k \ 
Remark 5. Theorem 4 is useful in multidimensional dynamical systems 
where involutivity of conservation laws can not be checked directly. 

7 Bi-Hamiltonian systems 

Further we will focus on non-Noether symmetries that satisfy condition (110). 
Besides yielding involutive families of conservation laws, such a symme- 
tries appear to be related to many known geometric structures such as bi- 
Hamiltonian systems and Frolicher-Nijenhuis operators (torsionless tangent 
valued differential 1-forms). The relationship between non-Noether symme- 
tries and bi-Hamiltonian structures was already implicitly outlined in the 
proof of Theorem 4. Now let us pay more attention to this issue. 
Originally bi-Hamiltonian structures were introduced by F. Magri in analisys 
of integrable infinite dimensional Hamiltonian systems such as Korteweg-de 
Vries (KdV) and modified Korteweg-de Vries (mKdV) hierarchies, Nonlin- 
ear Schrodinger equation and Harry Dym equation. Since that time bi- 
Hahmiltonian formalism is effectively used in construction of involutive fam- 
ilies of conservation laws in integrable models 

Generic bi-Hamiltonian structure on 2n dimensional manifold consists out of 
two Poisson bivector fields W and W satisfying certain compatibility condi- 
tion [W, W] = 0. If, in addition, one of these bivector fields is nondegenerate 
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(W n 7^ 0) then bi-Hamiltonian system is called regular. Further we will 
discuss only regular bi-Hamiltonian systems. Note that each Poisson bivec- 
tor field by definition satisfies condition (15). So we actually impose four 
restrictions on bivector fields W and W 



[W, W] = [W, W] = [W, W]=0 (126) 

and 

W n ^ (127) 

During the proof of Theorem 4 we already showed that bivector fields W and 
W = [E, W] satisfy conditions (126) (see (112)-(116)), thus we can formulate 
the following statement 

Theorem 5. Let (M, h) be regular Hamiltonian system on the 2n-dimensional 
manifold M endowed with regular Poisson bivector field W. Then, if the vec- 
tor field E on M generates the non-Noether symmetry, and satisfies condition 

[[E[E, W]]W] = 0, (128) 

the following bivector fields on M 

W,W — [E, W] (129) 

form invariant bi-Hamiltonian system ([W, W] = [W, W] = [W, W\ =0). 
Proof: See proof of Theorem 4. 

Sample. One can check that the non-Noether symmetry (53) satisfies con- 
dition (110) while bivector fields 

"-■k»k+frk (130) 



and 



d d d d 

W — [E, W] — Z\— — A — V z 2 — A — + 

OZ\ oz 3 oz 2 oz^ 



d d d d 

e23 " 24 7r- A 7T- + 7r- A 7r- ( 131 ) 

OZ 1 OZ 2 OZ 3 OZ4 
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form bi-Hamiltonian system [W, W] = [W, W\ = [W, W] = 0. Similarly, one 
can recover bi-Hamiltonian system of three particle Toda chain associated 
with symmetry (61). It is formed by bivector fields 




(132) 



and 




dz 3 dz 4 dz 4 dz 5 dz 5 dz 6 



(133) 



In terms of differential forms bi-Hamiltonian structure is formed by couple of 
closed differential 2-forms: symplectic form u> (such that dui — and uj n ^ 0) 
and oj* = Leuj (clearly duo' = dL^oJ = Lsdou = 0). It is important that by 
taking Lie derivative of Hamilton's equation 



along the generator E of symmetry 

L E (i Xh uJ + dh) = i[E,x h ]UJ + ix h L E uJ + L E dh = % Xh of + dL E h = (135) 
one obtains another Hamilton's equation 



where h* = L E h. This is actually second Hamiltonian realization of equa- 
tions of motion and thus under certain conditions existence of non-Noether 
symmetry gives rise to additional presymplectic structure u° and additional 
Hamiltonian realization of the dynamical system. In many integrable mod- 
els admitting bi-Hamiltonian realization (including Toda chain, Korteweg- 
de Vries hierarchy, Nonlinear Schrodinger equation, Broer-Kaup system and 
Benney system) non-Noether symmetries that are responsible for existence 



ix h u + dh = 



(134) 



i Xh uj' + dh' = 



(136) 
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of bi-Hamiltonian structures has been found and motivated further inves- 
tigation of relationship between symmetries and bi-Hamiltonian structures. 
Namely it seems to be interesting to know whether in general case exis- 
tence of bi-Hamiltonian structure is related to non-Noether symmetry. Let 
us consider more general case and suppose that we have couple of differential 
2-forms u> and u>* such that 

duj = du = 0, u n ^ (137) 



i Xh u) + dh = (138) 

and 

i Xh uf + dh' = (139) 

The question is whether there exists vector field E (generator of non-Noether 
symmetry) such that [E, = and ui* = L E oo. 

The answer depends on oA Namely if ui* is exact form (there exists 1-form 
9* such that oo* = d9') then one can argue that such a vector field exists and 
thus any exact bi-Hamiltonian structure is related to hidden non-Noether 
symmetry. To outline proof of this statement let us introduce vector field E* 
defined by 

i E .uj = 9' (140) 

(such a vector field always exist because uo is nondegenerate 2-form). By 
construction 

L E .u = lo' (141) 

Indeed 

L E »u! = di E 'U> + i E »duj = d9' = uj* (142) 

And 

i[E;x h ]V = L E .(i Xh uj) - i Xh L E .u = -d(E'(h) - h') = -dti (143) 
In other words [X h , E*\ is Hamiltonian vector field 

[X h ,E\=X h , (144) 
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One can also construct locally Hamiltonian vector field X g , that satisfies the 
same commutation relation. Namely let us define function (in general case 
it can be done only locally) 

g(z) = f h'dt (145) 
J 

where integration along solution of Hamilton's equation, with fixed origin 
and end point in z(t) = z, is assumed. And then it is easy to verify that 
locally Hamiltonian vector field associated with g(z), by construction, satis- 
fies the same commutation relations as E* (namely [X^,X g ] = Xy)- Using 
E* and Xy one can construct generator of non-Noether symmetry — non- 
Hamiltonian vector field E = E* — X g commuting with Xh and satisfying 

Lelu = Le»lu — Lx g uJ = Le*uj = ui* (146) 

(thanks to Liouville's theorem Lx g uJ = 0). So in case of regular Hamil- 
tonian system every exact bi-Hamiltonian structure is naturally associated 
with some (non-Noether) symmetry of space of solutions. In case when bi- 
Hamiltonian structure is not exact (u>* is closed but not exact) then due 
to 

to' = Leuj = (Heoj + ie^oj = d%E^> (147) 

it is clear that such a bi-Hamiltonian system is not related to symmetry. 
However in all known cases bi-Hamiltonian structures seem to be exact. 



8 Bidifferential calculi 

Another important concept that is often used in theory of dynamical sys- 
tems and may be related to the non-Noether symmetry is the bidifferential 
calculus (bicomplex approach). Recently A. Dimakis and F. Miiller-Hoissen 
applied bidifferential calculi to the wide range of integrable models includ- 
ing KdV hierarchy, KP equation, self-dual Yang-Mills equation, Sine-Gordon 
equation, Toda models, non-linear Schrodinger and Liouville equations. It 
turns out that these models can be effectively described and analyzed using 
the bidifferential calculi [17] [24]. Here we would like to show that each gen- 
erator of non-Noether symmetry satisfying condition [[-Ef-E", W"]] W] = gives 
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rise to certain bidifferential calculus. 

Before we proceed let us specify what kind of bidifferential calculi we plan 
to consider. Under the bidifferential calculus we mean the graded algebra of 
differential forms over the phase space 

n = U™ =0 Q (k) (148) 

(Vl^ denotes the space of /c-degree differential forms) equipped with a couple 
of differential operators 

d,d: n (k) -> n (k+1) (149) 

satisfying d 2 = d 2 = dd + dd = conditions (see [24]). In other words we 
have two De Rham complexes M,Q,d and M, Q, d on algebra of differential 
forms over the phase space. And these complexes satisfy certain compatibility 
condition — their differentials anticommute with each other dd+dd = 0. Now 
let us focus on non-Noether symmetries. It is interesting that if generator of 
the non-Noether symmetry satisfies equation W"]] W] = then we are 

able to construct an invariant bidifferential calculus of a certain type. This 
construction is summarized in the following theorem: 

Theorem 6. Let (M, h) be regular Hamiltonian system on the Poisson 
manifold M. Then, if the vector field E on M generates the non-Noether 
symmetry and satisfies the equation 

[[E[E,W]]W] = 0, (150) 

the differential operators 

du = $ u ([W,$ w (u)]) (151) 



du = $ u ([[E,W]$ w (u)]) (152) 

form invariant bidifferential calculus (d 2 = d 2 = dd + dd = 0) over the graded 
algebra of differential forms on M. 
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Proof: First of all we have to show that d and d are really differential 
operators , i.e., they are linear maps from f^ fc ) into satisfy derivation 

property and are nilpotent (d 2 = d 2 = 0). Linearity is obvious and follows 
from the linearity of the Schouten bracket [, ] and <&w, maps. Then, if u is 
a A;-degree form $vk maps it on /c-degree multivector field and the Schouten 
brackets [W, <liv(' u )] and [[E, W]$w( u )] result the k + 1-degree multivector 
fields that are mapped on k + 1-degree differential forms by $ w . So, d and d 
are linear maps from into f2( fc+1 \ Derivation property follows from the 
same feature of the Schouten bracket [, ] and linearity of <3>vk and $ w maps. 
Now we have to prove the nilpotency of d and d. Let us consider d 2 u 

d 2 u = 0> U ([W, $w($u,([W, *w(u)]))]) = * U ([W[W, ®w(u)}}) = (153) 

as a result of the property (112) and the Jacoby identity for [,] bracket. In 
the same manner 

dPu = ® w ([[W,E][[W,E}<S> w (u)}}) = (154) 

according to the property (116) of [W, E\ — W and the Jacoby identity. Thus, 
we have proved that d and d are differential operators (in fact d is ordinary 
exterior differential and the expression (151) is its well known representation 
in terms of Poisson bivector field). It remains to show that the compatibility 
condition dd + dd = is fulfilled. Using definitions of d, d and the Jacoby 
identity we get 

(dd + dd){u) = $ U ([[[W, E\W\<$> w {u)\) = (155) 

as far as (114) is satisfied. So, d and d form the bidifferential calculus over 
the graded algebra of differential forms. It is also clear that the bidifferential 
calculus d, d is invariant, since both d and d commute with time evolution 
operator W(h) = {h, }. 

Remark 6. Conservation laws that are associated with the bidifferential 
calculus (151) (152) and form Lenard scheme (see [24]): 

(k + 1) dl {k) = kdl (k+1) (156) 

coincide with the sequence of integrals of motion (98). Proof of this corre- 
spondence lays outside the scope of present manuscript, but can be done in 
the manner similar to [17]. 
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Sample. The symmetry (53) endows R A with bicomplex structure d, d where 
d is ordinary exterior derivative while d is defined by 



dz\ = z\dz\ — e z ' 3 ~ Z4 dz A 
dz 2 = z 2 dz 2 + e Z3 ~ Z4 dz 3 
dz 3 = z\dz 3 + dz 2 

dz A = z 2 dz A — dz\ (157) 

and is extended to whole De Rham complex by linearity, derivation property 
and compatibility property dd + dd = 0. By direct calculations one can 
verify that calculus constructed in this way is consistent and satisfies d 2 = 
property. To illustrate technique let us explicitly check that d 2 z\ = 0. Indeed 

d 2 zi = ddz\ = d(zidzi — e z ' A ^ ZA dz A ) = dz 1 A dz\ + z\ddz\ 
-e Z3 ~ Z4 dz 3 A dz A + e z ' s ~ Z4 dz A A dz A - e Z3 ~ Z4 ddz 4 = 
dz\ A dz\ — z\ddz\ — e Z3 ^ Z4 dz 3 A dz 4 

+e Z3 ~ Z4 dz 4 A dz A + e Z3 ~ Z4 ddz A = (158) 

Because of properties 

dzi A dz x = e Z3 ~ Z4 dzi A dz 4 , (159) 
-z x ddz x = z ie Z3 " Z4 dz 3 A dz A , (160) 
-e Z[, - Z4 dz 3 A dz A = -z l e Z3 ~ Zi dz l A dz A - e Z3 - Z4 dz 2 A dz A , (161) 
e Z3 ~ Z4 dz A A dz A = e Z3 ~ Z4 dz 2 A dz A (162) 

and 

e Z3 ~ Z4 ddz A = -e Z3 - Z4 d Zl A dz A (163) 
Similarly one can show that 

d 2 z 2 = d 2 z 3 = d 2 z A = (164) 
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and thus d is nilpotent operator d — 0. Note also that conservation laws 

J« = Zl + z 2 
/( 2 ) = z \ + z\ + 2e z *~ Z4 (165) 

form the simplest Lenard scheme 

2j/« = dl^ (166) 

Similarly one can construct bidifferential calculus associated with non-Noether 
symmetry (61) of three particle Toda chain. In this case d can be defined by 

dz\ = z\dz\ — e Z4 ~ Z5 dz$ 
dz 2 = z 2 dz 2 + e Z4 ~ Zr °dz A - e Z5 ~ Z(i dz 6 
dz 3 = z 3 dz 3 + e Z5 ~ Z6 dz 5 
dz^ = Zidz^ — dz 2 — dz 3 
dz$ = z 2 dz$ + dz\ — dz 3 

dz e = z 3 dz e + dz\ + dz 2 (167) 

and as in case of two particle Toda it can be extended to whole De Rham 
complex by linearity, derivation property and compatibility property dd + 
dd = 0. One can check that conservation laws of Toda chain 

=z 1 + z 2 

/( 2 ) = z \ + z \ + z \ + 2e Z4 ~ Z5 + 2e Zr '- Z6 
/( 3 ) = z \ + z\ + z\ + 3(zi + z 2 )e Z4 - Z5 + 3(z 2 + 2 3 )e 25 " 26 (168) 



form Lenard scheme 



2dl w = dl® (169) 
3d/ (2) = 2dl (3) (170) 
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9 Frolicher-Nijenhuis geometry 



Finally we would like to reveal some features of the operator R E (89) and 
to show how Frolicher-Nijenhuis geometry arises in Hamiltonian system that 
possesses certain non-Noether symmetry. From the geometric properties of 
the tangent valued forms we know that the traces of powers of a linear oper- 
ator F on tangent bundle are in involution whenever its Frolicher-Nijenhuis 
torsion T(F) vanishes, i. e. whenever for arbitrary vector fields X, Y the 
condition 



T(F)(X, Y) = [FX, FY] - F([FX, Y] + [X, FY] - F[X, Y]) = (171) 



is satisfied. Torsionless forms are also called Frolicher-Nijenhuis operators 
and are widely used in theory of integrable models, where they play role of 
recursion operators and are used in construction of involutive family of con- 
servation laws. We would like to show that each generator of non-Noether 
symmetry satisfying equation [[£'[£', W]]W] = canonically leads to invari- 
ant Frolicher-Nijenhuis operator on tangent bundle over the phase space. 
This operator can be expressed in terms of generator of symmetry and iso- 
morphism defined by Poisson bivector field. Strictly speaking we have the 
following theorem. 

Theorem 7. Let (M, h) be regular Hamiltonian system on the Poisson man- 
ifold M. If the vector field E on M generates the non-Noether symmetry 
and satisfies the equation 



then the linear operator, defined for every vector field X by equation 



is invariant Frolicher-Nijenhuis operator on M. 

Proof. Invariance of Re follows from the invariance of the Re defined by (89) 
(note that for arbitrary 1-form vector field u and vector field X contraction 
ixu has the property ir e xu = ixREU, so Re is actually transposed to Re)- 
It remains to show that the condition (110) ensures vanishing of the Frolicher- 
Nijenhuis torsion T(Re) of R E , i.e. for arbitrary vector fields X, Y we must 
get 



p[£,W]]W] = 



(172) 



R E (X) = ^ W (L E ^(X))-[E,X] 



(173) 
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T{R E )(X,Y) = [R E (X),R E {Y)\ - R E ([R E (X),Y] + [X,R E (Y)\ - 

R E ([X,Y]))=0 (174) 

First let us introduce the following auxiliary 2-forms 

to = $ W (W), to* = R E u uj" = ~R E uj* (175) 

Using the realization (151) of the differential d and the property (15) yields 

duo = $ U ([W,W]) = (176) 

Similarly, using the property (114) we obtain 

duo* = d<S> w ([E, W}) - dL E uj = & U ([[E, W]W}) - L E du = (177) 

And finally, taking into account that to' = 2& U) ([E, W\) and using the con- 
dition (110), we get 

du" = 2$ U ([[E[E, W]]W}) - 2dL E u' = -2L E du' = (178) 

So the differential forms lo,lo',lo" are closed 

duj = du' = duo" = (179) 
Now let us consider the contraction of T(R E ) and to. 

iT(R E )(X,Y) UJ = i[R E X,R E Y]U ~ i[R E X,Y]W* ~ i[X,R E Y]^' + i[X,Y] U " = 

L R E X^Y^* — iR E YLx^° — 

LR E xiYU* + i Y L RE xuj' - L x % ReY u* + i ReY L x u* + i[ X ,Y]v" = 

i Y L x u" - L x i Y u" + i [x , Y ]Cu" = (180) 

where we used (175) (179), the property of the Lie derivative 

L x i Y u = i Y L x u + i[x,Y]W (181) 
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and the relations of the following type 

Lr eX u = di RE xU + iR E xdu = di x uj' = L x uj' - i x du' = L x uj' (182) 

So we proved that for arbitrary vector fields X, Y the contraction of T(Re) (X, Y) 
and cu vanishes. But since W bivector is non-degenerate (W n ^ 0), its 
counter image 

to = $ U (W) (183) 

is also non-degenerate and vanishing of the contraction (180) implies that 
the torsion T(R E ) itself is zero. So we get 

T(R E )(X,Y) = [R E (X),R E (Y)} - R E ([R E (X), Y] + [X,R E (Y)\ - 

R E ([X,Y]))=0 (184) 



Sample. Note that operator R E associated with non-Noether symmetry 
(53) reproduces well known Frolicher-Nijenhuis operator 

d , d 

R E = z\dz\ ® azi <S> t; h 

az\ 0Z4 

9 7 d 

z 2 dz 2 <S> 7; h dz 2 <g> 7; h 

oz 2 oz 3 

d d 

zidz 3 ® — + e Zs Z4 dz 3 <g> — + 

OZ 3 0Z 2 

z 2 dz A <g> e Z3 ~ Z4 dz A <g> — (185) 

OZ A OZi 

(compare with [30]). Note that operator R E plays the role of recursion op- 
erator for conservation laws 

J« = z x + z 2 
/( 2 ) = z \ + z \ + 2e Z3 ^ 4 (186) 

Indeed one can check that 

2R E (dI {l) ) = dl {2) (187) 
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Similarly using non-Noether symmetry (61) one can construct recursion op- 
erator of three particle Toda chain 

R E = z x dz x <g> e Zi ~ Z5 dz 5 <g> — 

OZi OZ\ 

+z 2 dz 2 ® 7T- + e Z4 ~ Z6 dz 4 <g> 7T- - e Z5 ~ Z6 dz 6 <g> 

oz 2 oz 2 oz 2 

+z 3 dz 3 (g)^- + e Z5 ~ Z6 dz 5 <g> T f- 

+zidz 4 <g> rfz 2 ® - dz 3 <S> tt- 

0Z4 0Z4 0Z4 

9 1 9 J 9 

+z 2 dz 5 <g> h azi <g> dz 3 <g> 7— 

02:5 02:5 

<9 (9 

+z 3 dz 6 ® — + dz Y ® — + dz 2 ® — (188) 

C/2g 0^6 OZq 

and as in case of two particle Toda chain, operator R E appears to be recursion 
operator for conservation laws 

I W = z 1 + z 2 
/( 2 ) = z \ + z l + zl + 2e Z4 ~ Z5 + 2e Z5 ~ Ze 
J® = z \ + z \ + zl + 3( Zl + z 2 )e Zi - Z5 + 3(z 2 + z 3 )e Z5 ~ Z(i (189) 

and fulfills the following recursion condition 

rf/ (3) = 3R E (dI^) = Q(R E ) 2 (dI {1) ) (190) 



10 One-parameter families of conservation laws 

One-parameter group of transformations g a defined by (28) naturally acts on 
algebra of integrals of motion. Namely for each conservation law 

^J = (191) 
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one can define one-parameter family of conserved quantities J(a) by applying 
group of transformations g a to J 

J( a ) = g a (J) = e ahE J = J + aL E J + \(aL E ) 2 J + ... (192) 

Property (29) ensures that J(a) is conserved for arbitrary values of parameter 
a 

|j(-) = | fc W= fc (|j) = (193) 

and thus each conservation law gives rise to whole family of conserved quan- 
tities that form orbit of group of transformations g a . 

Such an orbit J(a) is called involutive if conservation laws that form it are 
in involution 

{J(a), J(b)} = (194) 

(for arbitrary values of parameters a, b). On In dimensional symplectic mani- 
fold each involutive family that contains n functionally independent integrals 
of motion naturally gives rise to integrable system (due to Liouville- Arnold 
theorem). So in order to identify those orbits that may be related to inte- 
grable models it is important to know how involutivity of family of conserved 
quantities J(a) is related to properties of initial conserved quantity J(0) = J 
and nature of generator E of group g a = e aLE . In other words we would like 
to know what condition must be satisfied by generator of symmetry E and 
integral of motion J to ensure that {J (a), J(b)} = 0. To address this issue 
and to describe class of vector fields that possess nontrivial involutive orbits 
we would like to propose the following theorem 

Theorem 8. Let M be Poisson manifold endowed with 1-form s such that 

[W[W(s),W](s)} = co[W(s)[W(s),W]] (c + -1) (195) 
Then each function J satisfying property 

W(L w(s) dJ) = Cl [W(s),W](dJ) ( Cl ^ 0) (196) 
(cq,i are some constants) gives rise to involutive set of functions 
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jM = {L w{s) ) m J { J {m \ J {k) } = 



(197) 



Proof. First let us inroduce linear operator R on bundle of multivector fields 
and define it for arbitrary multivector field V by condition 



R(V) = \([W(s), V] - <f> w (L w{8) $ u (V))) (198) 

Proof of linearity of this operator is identical to proof given for (89) so we 
will skip it. Further it is clear that 



R(W) = [W(s),W] (199) 

and 

R 2 (W) = R([W(s),W\) = l([W(s)[W(s),W}\ - <$> w {(L w{s) ) 2 u)) = 

^[WisWislW]} (200) 

where we used property 

$w((Lw(s)) 2 uj) = &w(Lw(s)Lw(s)u) = 
®w(iw(s)dL w{s) u}) + $ w (di W ( a )L W ( a )u) = [W^ w (i w{s) L w{s) u)} 

= [W[W(s),W}{s)}=c [W(s)[W{s),W]} (201) 

In the same time by taking Lie derivative of (199) along the vector field W(s) 
one gets 

[W[W(s), W](s)} = (L W{S) R + R 2 )(W) (202) 
comparing (200) and (202) yields 

(1 + c )(L w(s) R + R 2 ) = 2R 2 (203) 
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and thus 



L w{s) R = l —^R 2 (204) 
Further let us rewrite condition (196) as follows 

W(L w(s) dJ) = Cl R(W)(dJ) (205) 
due to linearity of operator R this condition can be extended to 

R rn (W)(L w{s) dJ) = Cl R m+1 {W){dJ) (206) 
Now assuming that the following condition is true 

W((L w(s) ) m dJ) = c m R m (W)(dJ) (207) 
let us take its Lie derivative along vector field W(s). We get 

R(W)((L w(s) ) m dJ) + W((L w(s) ) m+1 dJ) = 
mc m \^R m+1 (W) (dJ) + c m R m+1 (W) (dJ) (208) 

1 + Co 

where we used properties (199) and (204). Note also that (207) together with 
linearity of operator R imply that 

R k W((L w{s) ) m dJ) = c m R h+m (W)(dJ) (209) 
and thus (208) reduces to 

W((L w{s) r +1 d.J) = c m+1 R m+l (W)(dJ) (210) 
where c m+ i is defined by 

c m+ i = mc n ] C ° (211) 

i + Co 
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So we proved that if assumtion (207) is valid for m then it is also valid 
for m + 1, we also know that for m — 1 it matches (205) and thus by 
induction we proved that condition (207) is valid for arbitrary m while c n 
can be determined by 

c rn = c (m - 1)! 

Now using (207) and (209) it is easy to show that functions {L W ( s )) m J are 
in involution. Indeed 

{(L w{s) ) m J, (L w{s) ) k J} = W(d(L w{s) ) m JAd(L w{s) ) k J) = 

W((L w{s) ) m dJ A (L w(s) ) k dJ) = c m c k W(dJ A dJ) = (213) 

So we proved functions (197) are in involution. 

Further we will focus on concrete integrable models such as Toda chain, 
Broer-Kaup system and Benney system and we will use this theorem to prove 
involutivity of families of conservation laws constructed using non-Noether 
symmetries. 



11 Toda Model 

To illustrate features of non-Noether symmetries we often refer to two and 
three particle non-periodic Toda systems. However it turns out that non- 
Noether symmetries are present in generic n-particle non-periodic Toda chains 
too, moreover they preserve basic features of symmetries (53), (61). In case 
of n-particle Toda model symmetry yields n functionally independent con- 
servation laws in involution, gives rise to bi-Hamiltonian structure of Toda 
hierarchy, reproduces Lax pair of Toda system, endows phase space with 
Frolicher-Nijenhuis operator and leads to invariant bidifferential calculus on 
algebra of differential forms over phase space of Toda system. 
First of all let us remind that Toda model is 2n dimensional Hamiltonian 
system that describes the motion of n particles on the line governed by the 
exponential interaction. Equations of motion of the non periodic n-particle 
Toda model are 



1 - c 

1 +Cn 



III _L 



(212) 
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d 

dt q * = ^ 

= e (i - ly^-i* - € (n - i)e" +1 (214) 

(e(£r) = — e(— A;) = 1 for any natural k and e(0) = 0) and can be rewritten in 
Hamiltonian form (24) with canonical Poisson bracket defined by 

n d d 

corresponding symplectic form 

n 

uj = Y,dpiA d Qi (216) 
i=l 

and Hamiltonian equal to 



i n n—1 

^ = ^E^ + Ee" +1 (217) 
z i=i i=i 

Note that in two and three particle case we used slightly different notations 

Zi = Pi 

z n+l = qi i = 1,2,(3); n = 2(3) (218) 

for local coordinates. The group of transformations g a generated by the 
vector field E will be symmetry of Toda chain if for each Pi,qi satisfying 
Toda equations (214) g a {pi) , g a (qi) also satisfy it. Substituting infinitesimal 
transformations 

9a{Pi) =Pi + aE(pi) + 0(a 2 ) 

9a{Pi) = q l + aE{ qi ) + 0(a 2 ) (219) 
into (214) and grouping first order terms gives rise to the conditions 

jE( qi ) = Efa) 

jE( Pi ) = e(i - l) e *-^-*(S( ft _ 1 ) - E(q t )) - 

e ( n _ i) e «-*+i (E( qi ) - E(q i+1 )) (220) 
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One can verify that the vector field defined by 



E(pi) = -pi + e(i-l)(n-i + 2)e m - 1 ~ qi - e(n - i)(n - i)e qi ~ qi+1 + 
l -{t{% - \){vi-x +p i )e qi - 1 - qi - e(n - i)( Pi + p m )e^ +1 

i— 1 j n 

= (w - i + -^^Pk + ~ Pk + 

1 k=l 1 k=i+l 

-(Pi + e(i - \)e q ^- qi + e(n - i)e q ^ +1 ) (221) 

satisfies (31) and generates symmetry of Toda chain. It appears that this 
symmetry is non-Noether since it does not preserve Poisson bracket struc- 
ture [E, W] and additionally one can check that Yang-Baxter equation 
[[£'[£', W 7 ]]!^ 7 ] = is satisfied. This symmetry may play important role in 
analysis of Toda model. First let us note that calculating L E uo leads to the 
following 2-form 



n— 1 



L E oo = J2 Pi d Pi f\d qi + J2 ^ Ql+1 dqi A dq i+1 + ^ dp; A dpj (222) 

i<j 



i=l 



i=l 



and together to and L e oj give rise to bi-Hamiltonian structure of Toda model 
(compare with [30]). Thus bi-Hamiltonian realization of Toda chain can 
be considered as manifestation of hidden symmetry. In fact non-Noether 
symmetry carries even more information about the bi-Hamiltonian structure 
and give rise to the invariant symplectic potential for the differential form 
L e uj, i.e. invariant 1-form 9* such that 



d0' = u 



(223) 



This 1-form can be constructed by taking contraction of generator E of non- 
Noether symmetry and symplectic form u 



6' = t E u = J2 



^L2 



-Pidqi + (n-i + \)pidp { 



+ 



n-l 



e qi ~ qi+1 [(n-i+ l)dq i+1 - (n - t)dq t ] + - ^Tfec^ - Pj d Pi ) (224) 

i=l ^ i<j 
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In terms of bivector fields these bi-Hamiltonian system is formed by 



n 8 8 

'^§% A % (225) 



and 



n 8 8 ™ _1 8 8 8 8 

w = w] -S***is+S an ^ + g % A 

The conservation laws (45) associated with the symmetry reproduce well 
known set of conservation laws of Toda chain. 



n 

/ (D = C -(l) =S p. 

i=l 

n n—1 

7 (2) = (C7 (1) )2 _ 2C (2) = £ p 2 + 2 £ 

i=l i=l 
n n—1 

J(3) = C (l) )3 _ 3C -(1) C (2) + 3C r(3) = + 3 + pi+Oe*"^ 1 

2 = 1 1=1 



J(4) = C (l))4 _ 4(C (1))2 C7 (2) + 2(C ,(2) )2 + 4C -(1) C (3) _ 4C (4) = 
n n—1 

E Pi + 4 E (Pi + 2 P<Pi+i + P 2 i+i)e qi ~ qi+1 + 



i=i i=i 

n—1 n—2 



2 ^2 e 2 ( q '~ qi+1 ') + 4 E e q '~ qi+2 

i=l i=l 

m—1 

Jim) = (_l}rn+l mC (m) + ^ (^fc+lj^-fc)^) 

fc=l 

The condition [[^[-E, W]]W] = satisfied by generator of the symmetry E 
ensures that the conservation laws are in involution i. e. {C {k \C {m ^} = 0. 
Thus the conservation laws as well as the bi-Hamiltonian structure of the non 
periodic Toda chain appear to be associated with non-Noether symmetry. 
Using formula (88) one can calculate Lax pair associated with symmetry 
(221). Lax matrix calculated in this way has the following non-zero entries 
(note that in case of n — 2 and n = 3 this formula yields matrices (102)-(105)) 
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Lk,k — L>n+k,n+k — Pk 

L n +k,k+i = -L n +k+i,k = e(n - k)e qk ~ qk + 1 

m,k=l,2,...,n (228) 
while non-zero entries of P matrix involved in Lax pair are 



Pk,n+k 1 

p n+k k = - e (k - l) e qk - 1 - qk - e (n - k)e qk ~ qk+1 
P n+ k,k + i = e(n - k)e qk ~ qk ^ 
P n+ k,k-i = e(k - \)e qk --- qk 

fc = l,2,...,n (229) 

This Lax pair constructed from generator of non-Noether symmetry exactly 
reproduces known Lax pair of Toda chain. 

Like two and three particle Toda chain, n-particle Toda model also admits 
invariant bidifferential calculus on algebra of differential forms over the phase 
space. This bidifferential calculus can be constructed using non-Noether 
symmetry (see (152)), it consists out of two differential operators d, d where 
d is ordinary exterior derivative while d can be defined by 

dqi = pidqi + ^ dpj - dpj 

i<j i>j 

d Pi = pidpi - e qi - qi+1 dq i+1 + e qi ~'- qi d qi (230) 

and is extended to whole De Rham complex by linearity, derivation property 
and compatibility property dd + dd = 0. By direct calculations one can 
verify that calculus constructed in this way is consistent and satisfies d 2 = 
property. One can also check that conservation laws (227) form Lenard 
scheme 

(k + l)dl {k) = kdl (k+1) (231) 
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Further let us focus on Frolicher-Nijenhuis geometry. Using formula (173) 
one can construct invariant Frolicher-Nijenhuis operator, out of generator of 
non-Noether symmetry of Toda chain. Operator constructed in this way has 
the form 



Rt 



i=i 



9 j 9 
dpi<S)T^- + dqi <g) tt- 
oqi dpi 



n-l 



+ J2 z qi ~ x ~ 



u Pi i<j 



n-l 

-E* 

dpi C 



qi-qi- 


n dq i+ i ( 


d 

OPi 


d 


— dp j <8 


d ' 


d qj 


IT 



(232) 



One can check that Frolicher-Nijenhuis torsion of this operator vanishes and 
it plays role of recursion operator for n-particle Toda chain in sense that 
conservation laws 1^ satisfy recursion relation 



(k + l)R E (dI ik) ) = kdl {k+1) (233) 

Thus non-Noether symmetry of Toda chain not only leads to n functionally 
independent conservation laws in involution, but also essentially enriches 
phase space geometry by endowing it with invariant Frolicher-Nijenhuis op- 
erator, bi-Hamiltonian system, bicomplex structure and Lax pair. 
Finally, in order to outline possible applications of Theorem 8 let us study 
action of non-Noether symmetry (221) on conserved quantities of Toda chain. 
Vector field E defined by (221) generates one-parameter group of transfor- 
mations (28) that maps arbitrary conserved quantity J to 

J{a) = J + aJ« + ^ + ^ J(3) + ... (234) 

where 

jM = {L E ) m J (235) 

In particular let us focus on family of conserved quantities obtained by action 
of g a = e aLE on total momenta of Toda chain 



i=l 
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(236) 



By direct calculations one can check that family J(a), that forms orbit of 
non-Noether symmetry generated by E, reproduces entire involutive family 
of integrals of motion (227). Namely 



i n n—1 

= L E jW = {L E fj = 

in q n—1 

z j=i z i=\ 

= L E JW = {L E fj = 

n n—1 

( V n J - 3 Vdr 



o n n—1 

~ A Y.Pi + 3 E fe 2 + 2ftft+i + P* + i)e« +1 + 
i=i i=i 

q n—1 n— 2 

_ ^ g2(<Ji-9i+l) _|_ 3 ^ e 9i-9i+2 

^ i=l i=l 



J( m ) = LeJ^-V = (L E ) m J (237) 

Involutivity of this set of conservation laws can be verified using Theorem 8. 
In particular one can notice that differential 1-form s defined by 

E = W(s) (238) 

(where E is generator of non-Noether symmetry (221)) satisfies condition 

[W[W(s), W](s)\ = 3[W(s)[W(s),W}} (239) 

while conservation law J defined by (236) has property 

W(L w(s) dJ) = -[W(s),W](dJ) (240) 

and thus according to Theorem 8 conservation laws (237) are in involution. 



12 Nonlinear Schrodinger equation 

Toda model provided good sample of finite dimensional integrable Hamilto- 
nian system that possesses non-Noether symmetry. However there are many 
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infinite dimensional integrable Hamiltonian systems and in this case in order 
to ensure integrability one should construct infinite number of conservation 
laws. Fortunately in several integrable models this task can be effectively 
simplified by identifying appropriate non-Noether symmetry. First let us 
consider well known infinite dimensional integrable Hamiltonian system - 
nonlinear Schrodinger equation (NSE) 



where ip is a smooth complex function of (t,x) G R 2 . On this stage we will 
not specify any boundary conditions and will just focus on symmetries of 
NSE. Supposing that the vector field E generates the symmetry of NSE one 
gets the following restriction 



(obtained by substituting infinitesimal transformation ip — > tp + aE(ip) + 
0(a 2 ) generated by E into NSE). It appears that NSE possesses nontrivial 
symmetry that is generated by the vector field 



(here (f> is defined by <p x = ijii/j). 

In order to construct conservation laws we also need to know Poisson bracket 
structure and it appears that invariant Poisson bivector field can be defined 
if -0 is subjected zero ip(t,—oo) = ip(t, +00) = boundary conditions. In 
terms of variational derivatives the explicit form of the Poisson bivector field 
is 



ipt = i(ij>xx + 2^ 2 ^) 



(241) 



E{iP) t = i[E(i>) xx + 2i; 2 E(<i{j) + ±WE($)\ 



(242) 



E(ip) = i(ip x + -i) xx + ixj) + xip 2 ip) - t(ip xxx + 6ipi)ip x ) 



(243) 




(244) 



while corresponding symplectic form obtained by inverting W is 




■+oo 




(245) 



—00 
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Now one can check that NSE can be rewritten in Hamiltonian form 

rit = {hM (246) 
with Poisson bracket { , } defined by W and 

- ifj x ifj x ) (247) 

-oo 

Knowing the symmetry of NSE that appears to be non-Noether ([E, W] ^ 0) 
one can construct bi-Hamiltonian structure and conservation laws. First let 
us calculate Lie derivative of symplectic form along the symmetry generator 

/+oo 
[Stp x A Sijj + tf;8(f) A Sijj + tp8(f) A 5tp]dx (248) 
-oo 

The couple of 2-forms u and Leoj exactly reproduces the bi-Hamiltonian 
structure of NSE proposed by Magri [55] . Note also that using non-Noether 
symmetry one can construct invariant symplectic potential 

/+oo \ 
-oo Z 

-(V>V - M x )]dx + tS[ / (^ - M)dx] (249) 

L J— oo 

The the conservation laws associated with non-Noether symmetry are well 
known conservation laws of NSE 

/+oo 
-oo 

/+oo _ 
-oo 

J(3) = (C (1) )3 _ 3C (1) C (2) + 3C r(3) = 2 f +Q °^ 2 _ ^ x )dx 

J — oo 
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7 M = (_i)«>+i m( 7H + (-l) k+1 l( m -Vc {k) (250) 



k=i 



The involutivity of the conservation laws of NSE {C( k \ C*"^} = is related 
to the fact that E satisfies Yang-Baxter equation [[^[E 1 , = 0. So 

non-Noether symmetry of NSE reproduces infinite sequence of functionally 
independent conservation laws in involution, endows the phase space with 
invariant bi-Hamiltonian structure and gives rise to the following Frolicher- 
Nijenhuis operator 



This Frolicher-Nijenhuis operator plays the role of recursion operator for the 
infinite sequence of conservation laws 



13 Korteweg-de Vries equation 

Now let us consider other important integrable models - Korteweg-de Vries 
equation (KdV) and modified Korteweg-de Vries equation (mKdV). Here 
symmetries are more complicated but generator of the symmetry still can 
be identified and used in construction of conservation laws. The KdV and 
mKdV equations have the following form 




(251) 



dl^ = R E (dlW) 



(252) 



ut + u xxx + uu x = 0[KdV] 



(253) 



and 
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Ut + u xxx — 6u u x = 0[mKdV] (254) 

(here u is smooth function of (t,x) G R 2 ). The generators of symmetries of 
KdV and mKdV should satisfy conditions 



E{u) t + E{u) xxx + u x E(u) + uE{u) x = 0[KdV] (255) 

and 

E(u) t + E(u) xxx - 12uu x E(u) - 6u 2 E{u) x = 0[mKdV] (256) 

(again this conditions are obtained by substituting infinitesimal transforma- 
tion u — > u + aE{u) + 0(a 2 ) into KdV and mKdV, respectively). 
Further we will focus on the symmetries generated by the following vector 
fields 



E(u) = ^u xx + i« 2 + ^u x v + ^{u xxx + uu x ) - 
j^(6u xxxxx + 20u x u xx + 10uu xxx + bu 2 u x )[KdV] (257) 

and 



3t 
~2 



3 x 
(u xxxxx - 10u 2 u xxx - A0uu x u xx - 10u 3 x + 30u 4 u x )[mKdV} (258) 



(here v and w are defined by v x = u and w x = u 2 ) 

To construct conservation laws we need to know Poisson bracket structure 
and again like in the case of NSE the Poisson bivector field is well defined 
when u is subjected to zero u(t, — oo) = u(t, +oo) = boundary conditions. 
For both KdV and mKdV the Poisson bivector field is 



W = / dx— A — (259) 

J-oo OU OV 
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with corresponding symplectic form 



/+oo 
dx5u A Sv (260) 
-oo 

leading to Hamiltonian realization of KdV and mKdV equations 

u t = {h,u} (261) 

with Hamiltonians 

/+oo iy3 
(u 2 x - -)dx[KdV] (262) 
-oo O 

and 

/+oo 
(u 2 x + u 4 )dx[mKdV] (263) 
-oo 

By taking Lie derivative of the symplectic form along the generators of the 
symmetries one gets another couple of symplectic forms 

f+OO 2 

L e uj = / dx{5u A 5u x + -udu A 5v) [KdV] (264) 

J —oo O 

/+oo 
dx(5u A 5u x - 2u5u A Sw) [mKdV] (265) 
-oo 

involved in bi- Hamiltonian realization of KdV/mKdV hierarchies and pro- 
posed by Magri [55]. 

The conservation laws associated with the symmetries reproduce infinite se- 
quence of conservation laws of KdV equation 



I^ = C^ = 2 - ( +C °udx 

O J — oo 
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/( 2) = c m _ 2C( 2) = 4 2da . 

/(3) = (cd)) 3 - 3C«C( 2 ) + 3C^ = I / + °°(^ - u 2 x )dx 

9 J —oo o 

J(4) = ^(1)^ _ 4(C7 (1))2 C7 (2) + 2(C (2) )2 + 4^(1)^(3) _ 4C (4) = 

64 /■+«> 5 4 5 2 2 
45 J-oo 36 3 

m— 1 

7 M = (_!)™+i mC W + ^ (_!)*+!/("*-*)£(*) (266) 

k=i 

and mKdV equation 

/(i) =C (i) = _ 4 r 00 ^ 

J— oo 

/ + 00 
(« 4 + n 2 )cfe 
-oo 

/+oo 
(2« 6 + 10« 2 « 2 + 

J(4) = (C ,(1) )4 _ 4 ( C (1))2 C (2) + 2(C ( 2)) 2 + 4C (D C m _ 4C (4) = 

256 /"+ 00 

— / (5n 8 + 70w 4 « 2 - 7m 4 + Uu 2 u 2 xx + u 2 xxx )dx 

5 J —oo 

j(m) = (_ 1 )m+l mC( (m) + ^ (.^M-l^m-*)^*) 

fc=l 

The involutivity of these conservation laws is well known and in terms of the 
symmetry generators it is ensured by conditions [[-E[-E", W]]W] = 0. Thus the 
conservation laws and bi-Hamiltonian structures of KdV and mKdV hierar- 
chies are related to the non-Noether symmetries of KdV and mKdV equa- 
tions. Moreover these symmetries seem to be responsible for the existence of 
well known bi-Hamiltonian structures of KdV and mKdV equations. 



14 Nonlinear water wave equations 

Among nonlinear partial differential equations that describe propagation 
of waves in shallow water there are many remarkable integrable systems. 
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We already discussed case of KdV and mKdV equations, that possess non- 
Noether symmetries leading to the infinite sequence of conservation laws and 
bi-Hamiltonian realization of these equations, now let us consider other im- 
portant water wave systems. It is reasonable to start with dispersive water 
wave system, since many other models can be obtained from it by reduction. 
Evolution of dispersive water wave system is governed by the following set 
of equations 

u t = u x w + uw x 
v t = uu x — v xx + 2v x w + 2vw x 

Wt = w xx — 2v x + 2ww x (268) 

Each symmetry of this system must satisfy linear equation 

E{u) t = {wE{u)) x + {uE{w)) x 
E(v) t = (uE(u)) x - E(v) xx + 2{wE{v)) x + 2{vE{w)) x 

E(w) t = E{w) xx - 2E(v) x + 2{wE{w)) x (269) 

obtained by substituting infinitesimal transformations 

u — > u + aE(u) + 0(a 2 ) 
v -> v + aE(v) + 0(a 2 ) 
w + aE(w) + 0{a 2 ) (270) 

into equations of motion (268) and grouping first order (in a) terms. One 
of the solutions of this equation yields the following symmetry of dispersive 
water wave system 

E{u) = uw + x(uw) x + 2t(uw 2 — 2uv + uw x ) x 
3 

E(v) = -u 2 + Avw - 3v x + x(uu x + 2(vw) x - v xx ) 

+2t(u 2 w — uu x — 3v 2 + 3vw 2 — 3v x w + v xx ) x 
E{w) = w 2 + 2w x — Av + x(2ww x + w xx — 2v x ) 

—2t{u 2 + 6vw — w 3 — 3ww x — w xx ) x (271) 
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and it is remarkable that this symmetry is local in sense that E(u) in point 
x depends only on u and its derivatives evaluated in the same point, (this is 
not the case in KdV, mKdV and NLS equations where symmetries are non 
local due to presence of non local fields like v defined by v x = u in KdV 
equation, w defined by w x = u 2 in mKdV and (f> defined by (f> x = i/ji/j in case 
of nonlinear Scrodinger equation). 

Before we proceed let us note that dispersive water wave system is actually 
infinite dimensional Hamiltonian dynamical system. Assuming that u, v and 
w fields are subjected to zero boundary conditions 

u(±oo) = u(±oo) = w(±oo) = (272) 

it is easy to verify that equations (268) can be represented in Hamiltonian 
form 

u t = {h,u} 
v t = {h,v} 

w t = {h,w} (273) 

with Hamiltonian equal to 

\ r+oo 

h = — \ (u 2 w + 2vw 2 - 2v x w - 2v 2 )dx (274) 

4 J —oo 

and Poisson bracket defined by the following Poisson bivector field 



r+oo 


"1 5 


' 5 ' 


5 a 


' 6 ' 




w= / 










2^ 


Su 


5v 

X 


Sw 


.T- 



Now using our symmetry that appears to be non-Noether, one can calculate 
second Poisson bivector field involved in the bi-Hamiltonian realization of 
dispersive water wave system 



r+oo 


5 


' 5 ' 




u— A 




I — oo 


ov 


Su 



+ v— A 
ov 



Sv 



+ 



5v 



A 



Sw 



+ w— A 
ov 



Sw 



+ 



W 
_S_ 

Sw 



[E,W] 

A — 

Sw 
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Note that W give rise to the second Hamiltonian realization of the model 

u t = {h',u}. 
v t = {h',v}. 
w t = {h',w}, 

(277) 

where 

h* = — / [u 2 + 2vw)dx (278) 

4 J — CO 

and {, }. is Poisson bracket defined by bivector field W. 

Now let us pay attention to conservation laws. By integrating third equation 

of dispersive water wave system (268) it is easy to show that 

r+oo 

J (0) = / wdx (279) 

J — oo 

is conservation law. Using non-Noether symmetry one can construct other 
conservation laws by taking Lie derivative of J*- -* along the generator of 
symmetry and in this way entire infinite sequence of conservation laws of 
dispersive water wave system can be reproduced 

r+oo 

j(°) = / w dx 

J —oo 

J(i) = Le J(°) = -2 f^vdx 

/+oo 
{u 2 + 2vw)dx 
-oo 

/+oo 
{u 2 w + 2vw 2 - 2v x w - 2v 2 )dx 
-co 

= L E JW = (L E ) 4 JM = 

/+oo 
(u 2 w 2 + u 2 w x — 2u 2 v — 6v 2 w + 2vw 3 — 3v x w 2 — 2v x w x )dx 
-co 

J {n) = L E J {n ~ l) = {L E ) n J {0) (280) 

Thus conservation laws and bi-Hamiltonian structure of dispersive water 
wave system can be constructed by means of non-Noether symmetry. 
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To prove involutivity of infinite sequence of conserved quantities (280) one 
can use Theorem 8. In particular one can check that 1-form s defined via 

E = W(s) (281) 

(E is generator of non-Noether symmetry (271)) satisfies condition 

[W[W(s), W](s)\ = 3[W(s)[W(s), W}} (282) 

while 

/+oo 
vdx (283) 
-oo 

has property 

W{L w{s) dJ) = -[W(s),W](dJ) (284) 

and as a consequence of Theorem 8 gives rise to involutive family of conserved 
quantities. 

Note that symmetry (271) can be used in many other partial differential 
equations that can be obtained by reduction from dispersive water wave 
system. In particular one can use it in dispersiveless water wave system, 
Broer-Kaup system, dispersiveless long wave system, Burger's equation etc. 
In case of dispersiveless water waves system 



u t = u x w + uw x 
v t = uu x + 2v x w + 2vw x 

w t = —2v x + 2ww x (285) 

symmetry (271) is reduced to 

E(u) = uw + x(uw) x + 2t(uw 2 — 2uv) x 

3 

E{v) = -u 2 + 4vw + x{uu x + 2{vw) x ) + 2t(u 2 w — 3v 2 + 3vw 2 ) x 

E(w) = w 2 - Av + x(2ww x - 2v x ) - 2t{u 2 + 6vw - w 3 ) x (286) 
and corresponding conservation laws (280) reduce to 
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/-+oo 

J {0) = / wdx 

J -co 

/+oo 
-co 

/+CO 
{u 2 + 2^)cfe 
-co 

/+CO 
(u 2 w + 2vw 2 -2v 2 )dx 
-oo 

= L E J® = (L E ) 4 J(°) = 

/+oo 
(u 2 w 2 - 2u 2 v - 6v 2 w + 2vw 3 )dx 
-oo 

J in) = L E J {n - X) = {L E ) n J i0) (287) 



Another important integrable model that can be obtained from dispersive 
water wave system is Broer-Kaup system 

1 

1H = -Vxx + V X W + VW X 

w t = ~w xx + v x + ww x (288) 

One can check that symmetry (271) of dispersive water wave system, after 
reduction, reproduces non-Noether symmetry of Broer-Kaup model 



E(v) = Avw + 3v x + x(2(vw) x + v xx ) 
+t(3v 2 + 3vw 2 + 3v x w + v xx ) x 
E{w) = w 2 — 2w x + Av + x(2ww x — w xx + 2v x ) 

+t(6vw + w 3 — 3ww x + w xx ) x (289) 

and gives rise to the infinite sequence of conservation laws of Broer-Kaup 
hierarchy 



wdx 

oo 
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/'H-OO 
vdx 
-oo 
oo 

vwdx 

oo 

oo 



J(2) = Le J(1) = ( Le )2 J(0) = 4 y 

/+oo 
(to 2 + ^ + V 2 )<fe 
-oo 

JW = L£ J(3) = (L £ ) 4 J(°) = 

/ + O0 
(6v 2 w + 2to 3 + 3v x w 2 - 2v x w x )dx 
-oo 

J in) = L E J {n ~ x) = (L E ) n J^ 



(290) 



And exactly like in the dispersive water wave system one can rewrite equa- 
tions of motion (288) in Hamiltonian form 



v t = {h,v} 
w t = {h,w} 



where Hamiltonian is 

I r+oo 

h = - (to 2 + v x w + v 2 )dx 

1 J —oo 

while Poisson bracket is defined by the Poisson bivector field 



(291) 



(292) 



W 



r+oo 


' 5 


' 5 ' 




l—oc 


Sv 


5w 


X- 



dx 



(293) 



And again, using symmetry (289) one can recover second Poisson bivector 
field involved in the bi-Hamiltonian realization of Broer-Kaup system by 
taking Lie derivative of (293) 



W = [E, W] 



r+oo 


5 

v— A 


5 




5 


A 


5 


5 

+ w— A 


' 5 ' 




5 


' 6 ' 


1 — oo 


ov 


5v 


X 


5v 


X 


Sw 


ov 

X 


5w 


X 


5w 


5w 



^294) 
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This bivector field give rise to the second Hamiltonian realization of the 
Broer-Kaup system 



v t = {h',v}. 
w t = {h*,w}, 

(295) 



with 



| r+oo 

h* = — \ vwdx (296) 

4 J— oc 

So the non-Noether symmetry of Broer-Kaup system yields infinite sequence 
of conservation laws of Broer-Kaup hierarchy and endows it with bi-Hamiltonian 
structure. 

By suppressing dispersive terms in Broer-Kaup system one reduces it to more 
simple integarble model — dispersiveless long wave system 



v t = v x w + vw x 

w t = v x + ww x (297) 
in this case symmetry (271) reduces to more simple non-Noether symmetry 



E(v) = Avw + 2x(vw) x + 3t(v 2 + vw 2 ) x 
E(w) = w 2 + 4v + 2x(ww x + v x ) + t(6vw + w 3 ) x (298) 

while the conservation laws of Broer-Kaup hierarchy reduce to sequence of 
conservation laws of dispersiveless long wave system 



J (0) = / wdx 



J(i) = L E J^ = 2 J 
= L E JW = (L E ) 2 J^ = 4 J' 
J(3) = L E ,p) = (L E ) 3 JW = 12 f +0 °(vw 2 + v 2 )d. 



oo 
■hoo 

vdx 

oo 
oo 

vwdx 

oo 

+oo 

X 
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/+OC 
(3v 2 w + vw 3 )dx 
-oo 

j(») = L E j( n -V = (L E ) n jW (299) 

In the same time bi-Hamitonian structure of Broer-Kaup hierarchy, after 
reduction gives rise to bi-Hamiltonian structure of dispersiveless long wave 
system 



/+oo 
-oo 



5 

V— A 

ov 



Sv 



W 



+ w— A 
ov 



/+oo 
-oo 



'5 


' 5 ' 










5v 


Sw 


X- 



dx 



W = [E, W] 



Sw 



5w 



A 



_6_ 

Sw 



dx 



(300) 



Among other reductions of dispersive water wave system one should probably 
mention Burger's equation 



w t = w xx + ww x (301) 

However Hamiltonian realization of this equation is unknown (for instance 
Poisson bivector field of dispersive water wave system (275) vanishes during 
reduction). 



15 Benney system 

Now let us consider another integrable system of nonlinear partial differential 
equations — Benney system. Time evolution of this dynamical system is 
governed by equations of motion 

u t = vv x + 2(uw) x 
v t = 2u x + (vw) x 

w t = 2v x + 2ww x (302) 
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To determine symmetries of the system one has to look for solutions of linear 
equation 



E(u) t = (vE(v)) x + 2{uE{w)) x + 2{wE{u)) x 
E(v) t = 2E(u) x + (vE(w)) x + (wE(v)) x 
E(w) t = 2E(v) x + 2(wE(w)) x 



(303) 



obtained by substituting infinitesimal transformations 



u — > u -\- aE(u) + 0(a 2 ) 
v -> v + aE(v) + 0(a 2 ) 
w + aE(w) + 0(a 2 ) 



(304) 



into equations (302) and grouping first order terms. In particular one can 
check that the vector field E defined by 

E{u) = 5uw + 2v 2 + x(2(uw) x + vv x ) + 2t(Auv + v 2 w + 3uw 2 ) x 
E{v) = vw + 6u + x((vw) x + 2u x ) + 2t(Auw + 3v 2 + vw 2 ) x 

E(w) = w 2 + 4v + 2x(ww x + v x ) + 2t(w 3 + 4vw + Au) x (305) 

satisfies equation (303) and therefore generates symmetry of Benney system. 

The fact that this symmetry is local simplifies further calculations. 

In the same time, it is known fact, that under zero boundary conditions 



it(±oo) = i>(±oo) = «j(±oo) = 



(306) 



Benney equations can be rewritten in Hamiltonian form 



u t = {h,u} 
v t = {h,v} 
w t = {h,w} 



(307) 



with Hamiltonian 




(308) 
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and Poisson bracket defined by the following Poisson bivector field 



J-oo 2 5v 

Using symmetry (305) that in fact is non-Noether one, we can reproduce sec- 
ond Poisson bivector field involved in the bi-Hamiltonian structure of Benney 
hierarchy (by taking Lie derivative of W along E) 



5_ 

5v 



— A 

5u 



5 




5w 





dx 



(309) 



/ + OC 



5 

u— A 
ou 



5u 



5 

+ v— A 
ou 



5v 



5 

+ w— A 
ou 



5w 



W = [E, W] = 
5 



5v 



5w 



dttilO) 



Poisson bracket defined by bivector field W gives rise to the second Hamil- 
tonian realization of Benney system 

u t = {h',u}, 
v t = {h',v}. 
w t = {h',w}. 



with new Hamiltonian 



1 r+oo 

h' = - (v 2 + 2uw)dx 

6 J-oo 



(311) 



(312) 



Thus symmetry (305) is closely related to bi-Hamiltonian realization of Ben- 
ney hierarchy. 

The same symmetry yields infinite sequence of conservation laws of Ben- 
ney system. Namely one can construct sequence of integrals of motion by 
applying non-Noether symmetry (305) to 



J(0) 



wdx 



(313) 



(the fact that is conserved can be verified by integrating third equation 
of Benney system) . The sequence looks like 
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oc 
+00 

udx 

00 



f+00 

= / w dx 

J —00 

/■+00 

J« = L S J(°) =2 I vdx 
J(2) = = (L E ) 2 J^ = 8 f 

J — < 

/-(-OC 
(v 2 + 2u«;)cfa; 
-00 

= L E J^ = (L E ) A J^ = 

/+00 
(2uw 2 + Auv + v 2 w)dx 
-00 

J(5) = L £ jW = (L E f J(°) = 

/+00 
(4w 2 + 8uuu; + 2mu 3 + 2w 3 + v 2 w 2 )dx 
-00 

J (n) = W^ 1 ) = (L s ) n J (0) (314) 

So conservation laws and bi-Hamiltonian structure of Benney hierarchy are 
closely related to its symmetry, that can play important role in analysis of 
Benney system and other models that can be obtained from it by reduction. 
Similarly to the case of nonlinear water wave system, proof of involutivity of 
these conservation laws is based on Theorem 8. Namely, 1-form s defined by 
means of non-Noether symmetry (305) as follows 

E = W(s) (315) 

has property 

[W[W(s),W}(s)\ = \[W(s)[W(s),W}} (316) 
while conservation law 

/+00 
udx (317) 
-00 

satisfies condition 

W(L w(s) dJ) = -[W(s),W](dJ) (318) 

and thus according to Theorem 8 produces involutive family of conserved 
quantities. 
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16 Conclusions 



The fact that many important integrable models, such as Korteweg-de Vries 
equation, nonlinear Schrodinger equation, Broer-Kaup system, Benney sys- 
tem and Toda chain, possess non-Noether symmetries that can be effectively 
used in analysis of these models, inclines us to think that non-Noether sym- 
metries can play essential role in theory of integrable systems and properties 
of this class of symmetries should be investigated further. The present review 
indicates that in many cases non-Noether symmetries lead to maximal invo- 
lutive families of functionally independent conserved quantities and in this 
way ensure integrability of dynamical system. To determine involutivity of 
conservation laws in cases when it can not be checked by direct computations 
(for instance in direct way one can not check involutivity in many generic 
n-dimensional models like Toda chain and infinite dimensional models like 
KdV hierarchy) we propose analog of Yang-Baxter equation, that being sat- 
isfied by generator of symmetry, ensures involutivity of family of conserved 
quantities associated with this symmetry. 

Another important feature of non-Noether symmetries is their relationship 
with several essential geometric concepts, emerging in theory of integrable 
systems, such as Frolicher-Nijenhuis operators, Lax pairs, bi-Hamiltonian 
structures and bicomplexes. From one hand this relationship enlarges possi- 
ble scope of applications of non-Noether symmetries in Hamiltonian dynam- 
ics and from another hand it indicates that existence of invariant Frolicher- 
Nijenhuis operators, bi-Hamiltonian structures and bicomplexes in many 
cases can be considered as manifestation of hidden symmetries of dynam- 
ical system. 
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